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Abstract 

We study the transverse-momentum distribution of hadrons produced in semi-inclusive 
deep- inelastic scattering (SIDIS). We consider cross sections for various combinations of 
polarizations of the initial lepton and nucleon or the produced hadron, for which we per- 
form the resummation of large double-logarithmic perturbative corrections arising at small 
transverse momentum. We present phenomenological results for the processes Ip — > IttX 
with longitudinally polarized leptons and protons. We discuss the impact of the pertur- 
bative resummation and of estimated non-perturbative contributions on the corresponding 
cross sections and their spin asymmetry. Our results should be relevant for ongoing studies 
in the COMPASS experiment at CERN, and for future experiments at the proposed eRHIC 
collider at BNL. 
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1 Introduction 



Our knowledge about the structure of hadrons has been vastly improved by experiments 
with polarized high-energy lepton beams scattering off polarized nucleon targets. Spin ob- 
servables in deeply-inelastic lepton-nucleon collisions allow us to extract the spin-dependent 
parton distributions of the nucleon. At the same time, they challenge our understanding 
of the reaction mechanism within QCD, and our ability to perform reliable theoretical 
calculations of the relevant cross sections for each process and kinematic region of interest. 

Of particular interest in this respect are observables in semi-inclusive deeply-inelastic 
scattering (SIDIS), Ip —>■ IhX, for which a hadron h is detected in the final state. Depend- 
ing on the type of hadron considered, various different aspects of nucleon structure may 
be probed. Two current lepton scattering experiments, HERMES at DESY and COM- 
PASS at CERN, employ this method extensively. Measurements of spin asymmetries for 
longitudinally polarized beam and target, integrated over all transverse momenta of the 
produced hadron, have served to allow conclusions about the helicity-dependent up, down, 
and strange quark and anti-quark distributions in the nucleon [1]. It has also been recog- 
nized that distributions in the hadron's transverse momentum can be of great interest, in 
particular when the nucleon target is transversely polarized [2]. The associated experimen- 
tal investigations by HERMES [3], COMPASS [4], the SMC [5], and CLAS [6] have been 
remarkably productive and have opened windows on novel QCD phenomena such as the 
Sivers [7] and CoUins [8] effects. It is hoped that experiments at a possibly forthcoming 
polarized electron- proton collider, eRHIC, would carry on and extend these studies [9]. 
Theoretically, the most interesting kinematic regime is characterized by large virtuality 
of the photon exchanged in the DIS process, and relatively small transverse momentum, 
Qt ^ Q- This regime also provides for the bulk of the events in experiment. It is precisely 
here, for example, that effects related to intrinsic transverse momenta of partons in the 
nucleon may become visible, potentially offering new insights into nucleon structure. At 
the same time, the theoretical analysis of hard-scattering in this regime is fairly involved, 
but well-understood. In particular, the emission of gluons from the DIS Born process 
7*g — > q also leads to non-vanishing transverse-momentum of the final-state hadron and 
needs to be taken into account appropriately. It is the goal of this paper to present state-of- 
the-art calculations for the transverse-momentum dependence of some SIDIS observables. 
For this study, we will focus entirely on the set of leading-twist double-spin reactions. 



(i) 


e+p- 


^e + n + X , 


(ii) 


e +p - 


^e + A + X, 


(iii) 


e^-p^ - 


-^e + A^ + X , 


(iv) 


e + p- 


^e + n + X , 


(v) 


e+p - 


^e + A + X . 



Here arrows to the right (upward arrows) denote longitudinal (transverse) polarization. 
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Needless to say that the final-state pion could be replaced by any hadron. The same is true 
for the A, as long as the observed hadron is spin-1/2 and its polarization can be detected 
experimentally. In our study, we will make use of several ingredients available in the 
literature. In an earlier publication [10], two of us presented results for the 2^3 partonic 
reactions Iq — > Iqg and Ig — > Iqq, for all polarizations of interest. In the perturbative 
expansion and using coUinear factorization, these processes are the first to yield a non- 
vanishing transverse momentum of the produced hadron. In terms of the strong coupling 
as they are of order 0{as)- We therefore refer to them as "leading order (LO)" processes for 
the hadron transverse-momentum distribution. They are expected to be adequate (at least 
qualitatively) for achieving a good theoretical description at large transverse momentum, 
qT ~ Q. In the unpolarized case, the complete next-to-leading (NLO) {0{al)) corrections 
to the ^T-distribution have been calculated [11, 12, 13] which will lead to an improvement 
of the theoretical calculation at large g^. ^ 

At low transverse momentum, qj^ <^ Q, fixed-order calculations are bound to fail. The 
reason for this is well understood: when qx —>■ 0, gluon emission is inhibited, and the 
cancellation of infra-red singularities between real and virtual diagrams in the perturbative 
series leaves behind logarithmic remainders of the form 

, In- (gVg^) 

-2 l^J 

qrp 

in the cross section da/dq^ at the kth order of perturbation theory, where m = 1, . . . , 2A; — 1. 
Ultimately, when qx <^ Q, «s will not be useful anymore as the expansion parameter 
in the perturbative series since the logarithms will compensate for the smallness of ag. 
Accordingly, in order to obtain a reliable estimate for the cross section, one has to sum 
up ("resum") the large logarithmic contributions to all orders in a^. Techniques for this 
resummation are well established, starting with pioneering work mostly on the Drell-Yan 
process in the late 1970's to mid 1980's [18, 19, 20. 21, 22]. The "Collins-Soper-Sterman" 
(CSS) formalism [22] has become the standard method for qj^ resummation. It is formulated 
in impact-parameter (b) space, which guarantees conservation of the soft-gluon transverse 
momenta. The formalism has also been applied to the unpolarized SIDIS cross section [23, 
24], and it was found [24] that data for qx distributions from the HERA ep collider [25, 26] 
are satisfactorily described. 

Resummations at small transverse momentum q-r have also been developed for spin 
observables. In Refs. [27, 28] the CSS formahsm was apphed to longitudinal and trans- 
verse double-spin asymmetries in the Drell-Yan process. Early resummation studies on 
gT-distributions in jet production in polarized DIS were performed in [29]. In Ref. [30, 31] 
leading-logarithmic (LL) resummation effects were investigated for spin asymmetries that 
involve transverse-momentum dependent distributions, in particular the Collins functions 
mentioned above. In Ref. [32], resummation formulas for polarized SIDIS were derived, 

^The study [10] may be viewed as an extension of previous work on the qT-integrated polarized SIDIS 
cross section, for which the LO process is j*q q. LO calculations were performed in [14] and NLO ones 
in [15, 16] (for initial work on the NNLO corrections to the unpolarized gr-integrated SIDIS cross section, 
see [17]). 
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based on a factorization theorem at low transverse momentum [33] . A main result of [32] 
is that the CSS evolution equation is the same for the spin-dependent cases as in the unpo- 
larized one. Knowing the expression for the unpolarized resummed cross section, and the 
complete polarized LO cross sections [10], it is then relatively straightforward to determine 
the resummed expressions for the polarized case. These will be provided in explicit form in 
this paper, for the cases hsted in (1). We shall present results corresponding to resummation 
to next-to-leading logarithmic (NLL) accuracy, which corresponds to resummation of the 
towers with 2k — 1,2k — 2,2k — 3 in (2). We shall also present numerical estimates for the 
reaction ep — > enX at COMPASS and eRHlC, in order to study the general features of the 
resummation and its impact on the cross sections and the spin asymmetry. We note that 
in the numerical evaluation one needs to specify a recipe for treating the integration over 
the impact parameter at very large b ~ 1/Aqcd, in order to avoid the Landau pole present 
in the resummed expression. This is closely related to non-perturbative effects generated 
by resummation [22, 34, 35, 36, 37, 38, 39]. We will use two different methods for treating 
the large 6-region. 

We stress that we will not address single-transverse spin phenomena in this paper like 
those related to the Sivers and Collins effects mentioned above. For these, the analysis 
of QCD radiative effects is rather more involved [31] than for the double-spin case we 
consider, in particular regarding the connection of the behavior at small and large transverse 
momentum qt [40]. Also, as we shall see below, the gx-differential cross section in general 
contains terms depending on the angle (f) between the hadron and lepton planes. We shall 
only consider the resummation of the 0-independent pieces in the cross section, which 
dominate at small qx- It is an interesting topic by itself to study the resummation of the 
terms that depend on (p [41]. 

The remainder of this paper is organized as follows: Section 2 provides all formulas 
needed for the gr-differential SIDIS cross section, at LO and for the NLL resummed case. 
In section 3, we present numerical estimates for ep — > enX at COMPASS and eRHIC, based 
on the resummation formulas. We discuss in particular the effect of the resummation on 
the spin asymmetry. We also study the impact of the treatment of the large-6 region and 
of non-perturbative corrections on the resummed cross sections and the spin asymmetry. 
We conclude in Sec. 4. In the Appendix, we list the complete 0{as) cross sections for the 
processes in (1), correcting some typos in [10]. 

2 The gr- differential SIDIS cross section 
2.1 Kinematics 

We are interested in the cross section for the process 

tik) + A{pA, Sa) ^ l{k') + B{pB, Sb)+X , (3) 

where Sa, Sb are the spin vectors for the initial nucleon and the produced final-state hadron, 
which can be longitudinal or transverse, as indicated in (1). The lepton can be unpolarized 



4 



or longitudinally polarized; transverse-spin effects for the lepton are suppressed by nii/Q 
because of chirality conservation at the lepton-photon vertex and hence are negligible. Prom 
now on, we will for definiteness take hadron A to be a proton and the lepton / to be an 
electron. We define five Lorentz invariants, denoted Sep,Xiij,Q'^, zj, and q^, to describe the 
process. The center of mass energy squared, Sep, for the initial electron and the proton is 

Sep = (pa + kf c^2pA-k, (4) 

ignoring masses. The conventional DIS variables are defined in terms of the virtual photon 
momentum q — k — k'&s 

X,, = , g2 = -q' = -{k - k' f , (5) 

and they may be determined experimentally by observing the scattered electron. For the 
final-state hadron B, we introduce the scaling variable 

Pa ■ Pb f„s 

= ■ (6) 

PA-q 

Finally, we define the "transverse" component of q, which is orthogonal to both pa and pb'- 

tf = ,'-^rt-^pS. (7) 
PA -Pb Pa- Pb 

qt is a space-like vector, and we denote its magnitude by 

?T = y^. (8) 

To completely specify the kinematics, we need to choose a reference frame. We shall work 
in the so-called hadron frame [42, 43], which is the Breit frame of the virtual photon and 
the initial proton: 

?'' = (0,0,0,-Q), (9) 



Further, in this frame the outgoing hadron B is taken to be in the xz plane: 

pj = ££gfi + ^,%.o.t-iV (11) 



As one can see, the transverse momentum of hadron B is in this frame given by Zfqr- This 
is true for any frame in which the 3-momenta of the virtual photon and the initial proton 
are collinear. By introducing the angle between the hadron plane and the lepton plane, 
the lepton momentum can be parameterized as 

A;'* = ^ (cosh'0,sinh'0cos0,sinh'0sin0, — 1) , (12) 
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coseB= o smeB= , ■ (15) 



and one finds 

coshV'^^^^-l. (13) 

For the case (iii) in (1) [ep^ — > eA^^X] when the initial proton is transversely polarized 
and the transverse polarization of an outgoing spin- 1/2 hadron is observed, we need to 
parameterize the spin vectors Sa± and Sb± to lie in the planes orthogonal to pa and pb, 
respectively. In the hadron frame, they can be written as 

Sa± = (0,cos$A,sin$A,0) , 

3b± = (0,cos6BCos$B,sin$B,-sin6scos$B) , (14) 

so that ^A,B are the azimuthal angles of Sa,b± around 'Pa,b as measured from the hadron 
plane and 0b is the polar angle oi pb measured with respect to pa- One finds: 

We note that the polarization state depends on the frame we choose. For example, a state 
that is transversely polarized in the hadron frame becomes a mixture of longitudinally 
polarized and transversely polarized states in the laboratory frame where the initial electron 
and proton are coUinear. With the above definitions, the cross section for (3) can be 
expressed in terms of Sep, Xbj, Q^, Zf, and in the hadron frame. Note that (j) is 
invariant under boosts in the g*-direction, so that is the same in the hadron frame and, for 
example, in the photon-proton center-of-mass frame. 

2.2 Structure of the lowest-order cross section 

The complete set of the 0{as) spin-dependent partonic cross sections for large-g^- hadron 
production in SIDIS have been derived in Rcf. [10]. They are obtained from the Feynman 
diagrams for the reactions Iq Iqg and Ig Iqq. For completeness, and for the reader's 
convenience, we list them in the Appendix. Here we summarize the main characteristics of 
the hadronic cross section. 

The cross section can be decomposed into several pieces with different ^-dependences: 

for processes (i) and (n) in (1), 



^ ao + cos(0)cri , (17) 



dQ'^ dxbj dzf dq^dcf) 
for (iv) and (v), and 
d^(j 

dQHx,,dldqld<P = ^^^(^^ - - 2'^)'^' + - - + - ' 
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for (iii). In the hadron frame, cxo and cxj in these decompositions behave as as In (Q^/ q^) / 
at small qt- This is the manifestation of the large logarithmic corrections described in the 
Introduction at this order. At yet higher orders, corrections as large as a| In^'^ (Q^/ Qt) / Qt 
in the cross section arise. In the following, we will study the NLL resummation of these 
large logarithmic corrections within the CSS formalism. We note that the other components 
of the cross sections in Eqs. (16)-(18) are less singular as — ^ 0: ai and aj behave as 
asln(Q'^/q^) /qx, and 02 and crj as In (Q^/^f,). These pieces too, however, receive large 
higher-order corrections at small qT [23] and would need to be resummed [41]. In this paper 
we will only focus on the resummation for ctq and (Tq. 



2.3 Asymptotic part of the lowest-order cross section 

We will now discuss the behavior of the cross sections at small qr in more detail. This is a 
first step in arriving at the resummed cross section. Following [23, 24], we write the cross 
section at 0{as) as 



LO 



^5^asymp 



dQ^dxbjdzfdq^dcf) dQ'^dxbjdzfdq^dcf) 



Y{Sep,Q,qT,Xbj,Zf) , 



(19) 



where the first term on the right-hand-side is the "asymptotic" part in ctq and containing 
the pieces that are singular as In (Q'^/q^) /q^ or l/q^ at small qr, while Y collects the 
remainder which is at most logarithmic as qx ^ 0- It is straightforward to derive the 
asymptotic part of the LO unpolarized SIDIS cross section from the expressions in the 
Appendix. One finds: 



ofS^asymp 



2Cl 
Qt 



E 



dQ'^dxbjdzfdq^dcj) SnxljS^pQ'^ 



(20) 



where aem is the electromagnetic coupling constant, the fractional quark charge, and 
Cf — 4/3. Ai is defined in (A. 2). In Eq. (20), fq{x,n) and fg{x,n) denote, respectively, 
the unpolarized quark and gluon distribution functions for the proton at scale /i, and 
Dq{z, n) and Dg[z, n) are quark and gluon fragmentation functions for the produced hadron. 
Furthermore, for the unpolarized cross section, the P^J^'^^^ and Pj°"*''^°'* are identical, and they 

are equal to the customary LO unpolarized splitting functions: Pj* 
where 



m,(0) _ pout,{0) _ p(0) 



1 + 3^,^ 
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(21) 



Here, Tr — 1/2, and the "+" -prescription in the first fine acts in an integral from x to 1 as 



Ix ^(l — y)^- -fx 



^ , f{y) - / I , fn\^ n \ 
dy h/(l)ln(l -x) , 

i--y 



(22) 



for any suitably regular function /. Finally, the symbol ® denotes convolutions of the forms 



1 dx 



nn,{0) ( ^ 

\x 



(23) 



For the other processes listed in (1), the asymptotic parts are obtained in the same 
manner from the cross sections given in Appendix A. They can be cast in the form (20) 
with the following replacements: 
(ii) ep — >• eKX: 



Dr. 



AD 



9' 



(24) 



where the A/g, Afg and ADq, ADg denote the helicity-dependent quark and gluon parton 
densities and fragmentation functions, respectively (for definitions, see [2]). Furthermore, 



r,m,(0) 
ij 



AM 



(0) 



~,out,{0) 



AM 



(0) 



tj -I 



with 



APi? {x) = Pi? (x) , AP(J) (x) = [2 



x\ 



qq \ ^ qq 

(iii) ep^ eA^X: 

Ai sinh^ ■0 cos($A - $s - 20) 



AP<;j)(^) 



(25) 



Tr[2x-1] . (26) 



fq^ Sfq , /fl ^ , 



0, 



pi",(0) > ^p(O) po«t,(0) rp(0) 

aa 1 ao aa 



qq 



qq 



qq 



(27) 



Here, the 5fq and 5Dq denote the transversity distribution and fragmentation functions, 
respectively. As is well-known, gluons do not contribute to transversity at leading twist [2]. 
The LO transversity splitting function is given by 



5P^^\x)^Cf 



2x 

(T^ 



+ 



-r-5{l-x) 



(28) 
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(iv) ep eh.X: 



Al - A , U A/, , fg Afg , 



(v) ep eAX: 



,m,(0) 



(0) 



out,{0) 



(0) 



(29) 



Ai -Ae , 



D 



pm,(0) 
ij 



p(0) 



, AD,, 

rjOWt,(0) 



AD 



9 ' 



(30) 



2.4 NLL resummed cross section 



Resummation of the small- logarithms is achieved in impact-parameter (6) space, where 
the large corrections exponentiate. In 6-space, the leading logarithms q;| In^*^"^ (Q^/^t) /Qt 
turn into a^^Xxi'^ibQ /ho), where ho = 2e~^'^ with 7^; being the Euler constant. Subleading 
logarithms are down by one or more powers of ln(6(5/fco)- In the following, we will present 
the resummation formulas for the processes (i)-(v) in (1) we are interested in, based on the 
CSS formalism [22] (see also [23, 24, 32]). We will make some standard choices [22] for the 
parameters and scales in that formalism, which render the resulting expressions as simple 
and transparent as possible. 

For the unpolarized resummed cross one has 



■A 



/ 



dxbjdzfdQHqld(t> 2xlSyj (27r)2 
where the 6-space expression for the resummed cross section W is given by 



(31) 



Wib, Q, x,j, Zf) = Yl [Mx,j, bo/b) ® Cj^ 



C^"*®D,(^/,&o/&) ,(32) 



the convolutions being defined as in Eqs. (23). As indicated, the scale in the parton distri- 
butions is given by = bo/b. Wc will return to this in the next subsection. 
The Sudakov form factor S{b, Q) in Eq. (32) reads 



s{b, Q) = -l 



dkl 



(33) 



where the functions A and B have perturbative expansions of the form 



k=l 



a. 



TT 



k=l 



TT 



(34) 
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For the resummation at NLL, one needs the coefficients Ai^2 and Bi, which read [34, 44]: 



Ai — Cp , A2 — - Cp 



1 18 6, 



-Nf 
9 ^ 



(35) 



The functions Cj*" and Cj° in Eq. (32) are also perturbative. Their expansions read 



Qf/°"*(x,a,(/x)) = f2Cf^'^''^''\x) ' "^^^^ 

fe=0 



TT 



(36) 



Here // is a renormalization scale of order Q. The LO coefficients are given by 



19 ~ gg ~ ■ 

( in/out, (l) 



(37) 



The first-order coefficients (j^^/"'"^'^'-) could be obtained by expanding Eq. (32) to 0{as), 
performing an inverse Fourier transform, and comparing to the fixed-order {0{as)) Qt 
distribution. The asymptotic part we have given in Eq. (20) alone is not sufficient for 
this because it does not contain the contributions oc S'^{qT) at 0{as)- Fortunately, the 
coefficients are known in the literature, and have a very transparent structure in the MS 
scheme [24, 43] (see also work on the coefficients in the related Drell-Yan case in [21, 34, 45]). 
One has: 



(38) 



= —Pr^iz) + Hz)PP''"{z) + 5,,Cs5il - z) 



) out.e 



(39) 



Here the Pj*"'^(a;), /-"j™ '^(x) are the terms oc e in the LO splitting functions in 0? = 4 — 2e 
dimensions, P^^"'^^^^'^^\x,e), and are defined through 



-jin/out,{0) 
ij 



x) + eP, 



in/out,£ 
ij ' 



X) 



(40) 



The usual splitting functions in 0? = 4 dimensions have been given above in Eqs. (21). One 

has [46] 



Pgg^^) 



-C, 



\ — X — -^{^ ~ 2^) 



2Trx{\-x), P;^{x) = -Cfx. (41) 



Furthermore, the coefficient Cs in Eqs. (38), (39) arises from hard virtual corrections. It is 
therefore only present for the case ij = qq, where Cs = —ICp/^- Finally, for the "final- 
state" coefficient C°^^'^^^ there is an extra contribution oc In^; [24, 43] times the LO splitting 
function. It is of kinematic origin and contributed by the phase space for coUinear gluon 
radiation off the quark leg in the final state [49] . 



10 



With this, the formulas for the resummation of the unpolarized SIDIS cross section 
in 6-space are complete. It is now straightforward to extend them to the various spin- 
dependent cross sections in (1). The Sudakov form factors are entirely related to soft-gluon 
emission and hence are spin-independent and the same in each case [32]. For each process, 
we first need to make the appropriate replacements as described above in Eqs. (24)- (30). In 
addition, we obtain the respective coefficients from Eqs. (38), (39), substituting 

appropriately all splitting functions, and also the p.J^/"'^^'^ by the cx e terms in the o? = 4— 2e- 
dimensional versions of the LO splitting functions that (24)-(30) direct us to use. For this we 
need the oc e terms in the spin-dependent splitting functions. For longitudinal polarization 
we have [47] 



AP;,(a;) = -Cf 



l-x--6(l 
2 ^ 



X] 



AP^lix) = -2Tn{l - x) , AP;,(x) = 2C^(1 - x) . 



91^ 



(42) 



In case of transversity, the splitting function in 4 — 2e dimensions is for a; < 1 identical to 
that in four dimensions [48], so that 



Cf 



6{l-x) 



(43) 



The coefficient for the case ij = qq in Eqs. (38), (39) is the same for all polarized cases, 
since it comes from virtual diagrams. 

Knowledge of the coefficients Ai^2, Bi and C\f'^^'' in Eqs. (31), (32) is sufficient for 
the resummation to NLL accuracy. For consistency, we also need to use the two-loop 
expression for the strong running coupling, and the NLO evolution for the parton densities 
and fragmentation functions. Let us finally give an explicit formula for the NLL expansion 
of the Sudakov form factor. Defining 



L = In (l + g , 

and choosing the renormalization scale we have [50]: 

1 



k{ois{ij)L) fi{as{ij)L) , 



where 



/o(y) 
Hv) 



A 

1 

■Kf3o 



[l3oy + \n{l-Poy)]: 



(44) 



(45) 



(46) 



■In^d 



PoV ^ ln(l - Poy) 

1- Pov 1 - Poy 



+ 



■Kpo 



ln(l - Poy) 




ln(l - Poy) + 



Poy 



i-Poy 



(47) 
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with the coefficients of the QCD /3-function Po = (33 - 2iV/)/(127r) and Pi = (153 - 
19A^F)/(247r2). We observe that as y ^ 0, fo{y) = 0{y'^) and fi{y) = 0{y). Note that in 
Eq. (44) we have followed Ref. [51] to choose In (1 + Q%'^/hl) rather than \n.{Q%'^ /hi) as 
the large logarithm in 6-space. As far as the large-6 behavior and NLL resummation are 
concerned, the two are equivalent, of course. With the choice in (44), however, the Sudakov 
exponent vanishes at 6 = 0, as it should [19, 21, 22], whereas for the other choice it gives 
large logarithms not only at large h (small qt) but also artificially at small h (large qt). For 
further discussion, see [51]. 

2.5 Evolution of part on distributions and fragmentation 
functions 

As we have seen in Eq. (32), in the CSS formahsm the parton distributions and fragmenta- 
tion functions are evaluated at the factorization scale = ho/h. As discussed in [37], it then 
becomes a great convenience to treat them in Mellin-moment space, because this enables 
one to explicitly express their evolution between a large scale ~ Q and ho/h. In this way, 
one avoids the problem normally faced in qr resummation that one needs to call the parton 
densities or fragmentation functions at scales far below their range of validity (see below), 
so that some sort of "freezing" (or related prescription) for handling them is required. We 
also anticipate that below we will choose the impact parameter h to be complex-valued, so 
that it becomes desirable to separate the complex scale ho/h from the scale at which the 
parton densities are explicitly evaluated. 

To be more specific, we write the resummed cross section as W{h, Q, x^j, Zf) in Eq. (32) 
as an inverse Mellin transform of moments in xi,j and Zf: 

W{b, Q, x,„ Zf) = E (^) ' dN j^^ dM zj^x^/' 

j = <1,<1 
hk = q,q,g 

X (6o/&)C^' V(^'«)C^"*'^L>f (6o/&) , (48) 
where the moment of each function is defined as 

(/.) ^ x''-'Mx,^l) , C;r^''{a,i^l)) ^ [ a;^-^C;T(x, «,(/.)) , (49) 

and so forth. The integrations over and M in Eq. (48) are over contours in the complex 
plane, to be chosen in such a way that they lie to the right of the rightmost singularities of 
the parton distributions and the fragmentation functions. In moment space, it is possible 
to express the moments //^(feo/&) and Djf(bo/b) by the corresponding moments at scale 
~ Q. The evolution "factor" between the scales bg/b and Q (which in general is a matrix 
because of singlet mixing) can be written in closed analytical form and can be expanded to 
NLL accuracy in L. One then only needs the parton densities and fragmentation functions 
at scale Q. For further details, see [37]. 
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2.6 Inverse Fourier transform and non-perturbative corrections 



As shown in Eq. (31), the rcsummcd cross section in gj--spacc is obtained by an inverse 
Fourier transform of W{b,Q,Xhj, Zj). This involves an integration over b = \b\ from to 
oo. Because of the Landau pole of the perturbative strong coupling in the Sudakov expo- 
nent (33) at kr = Aqcd, this integration is ill-defined. In the expansion (45) for the Sudakov 
form factor, the singularity occurs at f3Qas{n)L = 1. In other words, the 6-integration ex- 
tends over both perturbative {b ^ 1/Aqcd) and non-perturbative {b < 1/Aqcd or even 
larger) regions. In order to define a perturbative resummed cross section, a prescription for 
the b integration is required that avoids the Landau pole. The ambiguity in the perturbative 
series implied by the presence of the Landau pole corresponds to a non-perturbative cor- 
rection. One may therefore use resummation to examine the structure of non-perturbative 
corrections in hadronic cross sections [52]. Typically, for the gy-differential Drell-Yan or 
SIDIS cross sections, one expects Gaussian non-perturbative effects in 6-space [22, 39]: 

g5(6,Q) _^ ^S{b,Q)-gb^ ^ ^^q^ 

where g is a. coefficient that may be determined by comparison to data. It will depend 
in general on Q and on Xbj and Zf. This non-perturbative Gaussian form factor may be 
thought of to partly incorporate (or, at least, mimic) the effect of smearing due to the 
partons' intrinsic transverse momenta. 

Several prescriptions for treating the inverse Fourier transform have been proposed in 
the literature. The earliest one, known as "6*-prescription" [22, 34], is to prevent b from 
becoming larger than a certain b^^x in the Sudakov form factor by replacing b ^ b* — 
b/^1 + &^/&max- This evidently introduces a new parameter, 6max- The prescription was 
refined recently in [36]. Extensive phenomenological studies using this prescription and a 
non-perturbative term as in (50) have been carried out in particular for the gj-distribution 
in the Drell-Yan process [34, 36, 35]. The method has also been used in studies of the 
gT-distribution in unpolarized SIDIS [24, 43, 53]. 

Another method, proposed in [54] and apphed to the Drell-Yan cross section in [37] , is 
to deform the b integration in Eq. (31) into a contour in the complex 6-plane. In this way, 
the Landau pole is avoided since it lies far out on the real-6 axis. This procedure does not 
introduce any new parameter and is identical to the original 6-integration in (31) for any 
finite-order expansion of the Sudakov exponent. For all details of the complex-6 prescription, 
see [37]. This method was also recently used in a study of the transverse-momentum 
distribution of Higgs bosons at the LHC [51, 55]. In the present paper we will use both 
the b* and the complex-6 methods and compare the results and their dependence on the 
non-perturbative parameters chosen. We note that for the complex-6 method the 6-integral 
converges even for g = in Eq. (50), while for the b* prescription a (non-perturbative) 
suppression of the integrand is required [37]. We finally mention that other methods for 
dealing with the behavior of the resummed exponent at large b have been discussed and 
used [38]. Also resummations directly in gr-space have been studied [56]. 
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2.7 Matching to finite order 



In order to obtain an adequate theoretical description also at large qt ^ Q, we "match" the 
resummed cross section to the fixed-order (LO, 0{as)) one. This is achieved by subtracting 
from the resummed expression in Eq. (31) its expansion to 0{as), 

and then adding the full 0{as) cross section, given by Eq. (19). We use this matching 
procedure, which avoids any double-counting of higher perturbative orders, in an identical 
way for the various polarized cross sections. 



3 Numerical analysis for ep enX 

In order to study the effect of resummation on the polarized SIDIS cross section and 
spin asymmetry, we will perform a numerical calculation for one case considered in (1), 

— * 

Ip — > IttX. This calculation is relevant for the ongoing COMPASS experiment and for 
experiments at the proposed polarized ep collider cRHIC. As typical values of the kincmat- 
ical parameters, we choose S^p = 10^ GcV^, = 100 GeV^, x^j = 0.012 for eRHIC and 
S^p = 300 GeV\ g2 _ iQ QgY2^ ^ q Q4 COMPASS. Clearly these choices are only 
meant to represent the typical kinematics; future detailed comparisons to experimental data 
will likely require to study a much broader range of values. We integrate over Zf > 0.2. 
Our two sets of the parameters give an identical cosh-?/' in Eqs. (13) and (A. 2). For the 
NLO parton densities and the fragmentation functions, we use the GRV unpolarized distri- 
butions [57], the GRSV "standard" polarized distributions [58], and the pion fragmentation 
functions by Kretzer [59] . We have each of these available in MeUin-moment space, which 
is very helpful in the light of the discussion in Subsec. 2.5. 

In addition, we will use a simple Gaussian non-perturbativc factor as in Eq. (50). We 
emphasize again the more illustrative character of our study and will therefore choose only 
some representative values for g, in order to investigate the sensitivity of the results to g. 
For the complex-6 prescription, we will use g — 0.6, 0.8 GeV^ for both the eRHIC and 
the COMPASS cases. The value g = 0.8 GeV^ was found in [37] in the context of the 
"joint" resummation for Z production at the Tevatron. This value therefore applies to the 
scale Mz- g is expected to have logarithmic dependence on Q [35] and should be smaller 
at lower energy. This motivates our lower choice of = 0.6 GeV^. For the b* method, 
we will choose bmax — l/('\/2 GeV) and the values g — 0.8, 1.3 GeV^ for the eRHIC case 
and g = 0.4, 0.8 GeV^ for COMPASS. We note that in the analysis [24, 43] of the HERA 
data [25, 26] for SIDIS observables, a detailed phenomcnological study of non-perturbative 
effects was performed, using the b* prescription. From the energy flow observable in SIDIS 
it was found that the non-perturbative effects appear to become larger at small x^j. If we 
boldly assumed that the non-perturbative form of [24, 43] also applied for the kinematics 
and the observable we are considering here, we would find g 0.4 GeV^ for the COMPASS 



14 




60 



'o 40 



20 







1 1 1 1 

(c) Spin asymmetry 




S=100^ GeV^ 




Q^=100 GeV^ 




=0.012 




z, -integrated (z, >0.2) 












- . _ . 




\ 


/ 


- - LO 




— - Complex-b (g=0) 


- 1 


Complex-b (g=0.6) 


1 


Complex-b (g=0.8) 




b*-metliod {g=0.8) 


1 

1 


b*-metliod (g=1.3) 

1 1 1 



2 4 6 8 10 

qj[GeV] 



Figure 1: (a) Unpolarized SIDIS cross section for eRHIC kinematics. We show the fixed- 
order (LO) result, and resummed results for the complex-6 method with non-perturbative 
parameters g = and g = 0.6, 0.8 GeV^, and for the b* method with 6max = l/('\/2 GeV) 
and g = 0.8, 1.3 GeV^. (b) Same for the longitudinally polarized case, (c) Spin asymmetries 
corresponding to the various cross sections shown in (a) and (b). 

situation, and g ^ 1.3 GeV^ at eRHIC. We are also interested in the comparison of the 
two methods with the same non-perturbative Gaussian, hence the additional choice of 
g — 0.8 GeV^. Finally, we shall always use the same values of g in the unpolarized and the 
polarized cases. This assumption is motivated by the fact that the perturbative Sudakov 
form factors are independent of polarization. 
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Figures 1 (a) and (b), respectively, show the unpolarized and polarized cross sections 



2^ LI ^'^^'^ 



d(A)cT 



(52) 



dxbjdzfdQ'^dqTd4> 



for eRHIC kinematics, where z^'^^ is given in Eq. (A. 6). Note that we have multiphed 
by a factor qx, as compared to the cross section we considered in (19). The LO cross 
sections rapidly diverge as Qt 0. For the matched cross section using the complex-6 
method with g = 0, one obtains an enhancement at intermediate qt and the expected 
reduction at small qt- The inclusion of the non-perturbative Gaussian form factor makes 
this tendency stronger. However, the results for different choices of the Gaussian, g = 0.6 
and 0.8 GeV^, are not very different and just have a slightly smaller normalization and are 
shifted to the right. Also shown in these figures are the curves for the b* prescription with 
&max = 1/(^2 GeV) and g = 0.8 and 1.3 GeV^. For the eRHIC case, the choice of ^ = 0.8 
gives a result relatively close to the one with = in the complex-6 method, while g — 1.3 
gives a result very close to those with g — 0.6 and 0.8 in the complex-6 method. One should 
note that these resummed curves all give a very similar g^-integrated cross section, close to 
the full NLO one, due to our choice of L in Eq. (44) and our matching procedure described 
in the previous subsection [51]. 

Figure 1 (c) shows the corresponding spin asymmetries, defined by the ratios of the 
polarized and unpolarized cross sections for all the curves shown in Figs. 1 (a) and (b). 
Although the effects of rcsummation and the non-perturbative Gaussians are significant in 
each cross section, they cancel to a large degree in the spin asymmetry. If one looks in more 
detail, resummation somewhat enhances the asymmetry compared to LO in the range of 
small to intermediate qr, where also the cross sections shown in Figs. 1 (a) and (b) receive 
enhancements due to the resummation. 

Figures 2 (a)-(c) show the same quantities as Figs.l (a)-(c), now for the COMPASS 
kinematics described above. In this case, the complex-6 method without Gaussian smearing 
{g = 0) turned out to be difficult to control numerically at small qt, and we do not show 
the result for it. We find that resummation leads to a significant enhancement of the cross 
section at g^ > 1 GeV. In both unpolarized and polarized cross sections, the resummed 
results we show, for the complex-6 method with g = 0.6 and 0.8 GeV^, and for the b* 
prescription with 6max = 1/(a/2 GeV) and g = 0.8 GeV^, turn out to be very similar, 
while the b* prescription with g = 0.4 GeV^ gives a higher peak that is shifted to the left 
compared to the other three resummed results. All these resummed results give a very 
similar spin asymmetry for the process Ip — > InX for COMPASS kinematics; we find that 
resummation just leads to a moderate decrease of the asymmetry. 



We have carried out a study of the soft-gluon resummation for the transverse-momentum 
{qt) distribution in semi-inclusive deeply-inelastic scattering. Resummation is crucial at 
small transverse momenta, qt <^ Q, where it takes into account large double-logarithmic 



4 Summary and Conclusions 
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(b) Polarized cross section 
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Figure 2: Same as in Figs. 1 (a)-(c), but for COMPASS kinematics. For the b* prescription, 
we have chosen here the non-perturbative parameters g — 0.4, 0.8 GeV^. 

corrections to all orders in the strong coupling constant. We have considered all relevant 
leading-twist double-spin cross sections, focusing on the terms that are independent of the 
angle between the lepton and the hadron planes, and have presented the resummation 
formulas for each. 

— * 

We have performed phenomenological studies for the process Ip InX at COMPASS 
and at a possible future polarized ep collider, eRHIC. Here we have chosen two different 
prescriptions for treating the region of very large impact parameters in the Sudakov form 
factor, which is related to the onset of non-perturbative phenomena. We have used simple 
estimates for the non-perturbative term suggested by the resummed formula. Our results 
indicate that resummation effects as well as non-perturbative effects cancel to a large extent 
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in the spin asymmetry. 
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A Analytic formulas for the LO cross sections 

Here we summarize the LO cross section formulas for the processes in (1), which were 
derived in [10]. The differential cross sections are given by 



rmn 



OC J Zm.i.T). ^ 



where the Ak are defined as [10, 42] 

^1 = 1-1- cosh^ , 
A2 = -2, 

^3 = — cos (f) sinh 2t/j , 
^4 = COS 20 sinh^ ip , 
Ae — —2 cosh ■0 , 
At = 2 cos 4> sinh ip , 
^8 = — sin (f) sinh 2ip , 

Ag = sin 20 sinh V, (A. 2) 

and the summation over /c is A; = 1, . . . , 4 for processes (i) and (ii) in (1), A; = 1, . . . , 4, 8, 9 
for process (iii), and A; = 6, 7 for processes (iv) and (v) (see below), aem — e^/An is the 
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QED coupling constant, and we have introduced the variables 



Xbj ^ Zf 

X — — - , z — — 
X z 



and 



Xr, 



(A.3) 



(A.4) 



For a given Sep, and qt, the kinematic constraints for Xbj and Zf are 



< Xbj < 1 , 



< Z < 1 - Xbj 

1 - Xbj + XbjqT/Q"^ 



(A.5) 
(A.6) 



Consequently, qt is limited by 



1 



\ \^bj ) \Zf 
In the hadron frame, the transverse momentum, pr, of p'^ obeys 



(A.7) 



PT = ZfQT < ZfQ, 



\ i^bj {zf ^) ■ 



(A.8) 



The term [/ o D o dk] in Eq. (A.l) takes the following form for the processes in (1): 
{i)e+p^e + n + X [60]: 

[foDoa,] = Y: e%{x)D,{z)a'^^ + Y: elf,{x)D,{z)dl + ^ , (A.9) 



9,9 



9,9 



9,9 



where {Cp = 4/3) 



— 2al — SCpXZ , 



99 



(A.IO) 



19 



^Ig = 2ctJ^ = 8^(1 - x) , 
(ii) e + p^ e + A + X: 



where 



Q,Q ?,? 

+ J2elAf,{x)AD,{z)A,a'^^, 



91 

1,1 



Al5J, = 5J, (/c = 1,2,3,4), 

, ^ {2x - 1) - ly - g^x^} 

^"^^^ Q^q^x{x - 1) 

A^aJ, = A^af = , 

.3 2{Q^(x-l)-g^£} 
Alct = — , 



X - 1 ' 



3 4CFx2j ..2 ,^^(1t\(1t 

20 



(iii) e + p^ ^ e + + X: 

For this process, it is more transparent to write the cross section by including the factors 
Ak in (A.l): 



fe=l, ■••,4,8,9 



(A.17) 



Q,Q 



where 



Q 



(l + cosh^ i/j) cos($A - - — sinh 2'0 cos($a - $b - 0) 



sinh ijj cos($^ - $s - 20) 
Qt 



(A.18) 



Here the terms with ^3 and^s in Eq. (A.l) have been combined to give the cos — 
term in (A.18). Likewise, those with ^4^9 give the term with cos($a — $b — 20). 
(iv) e + p — > e + TT + X: 



where 



[foDoa,]=Y: elAf,ix)D,{z)ALoB',, + E e^A/,(x)r',(z) A^o^J, 

+ EejA/,(a;)D,(;^)A^o5j, , (A.19) 



9,9 



2(:7p<'(^ + S2)?- 



xz 



Qt 



xzq^ 1 



(A.20) 
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2x - 1 / ^ X - 1 g2 



2Q (x- 1)(2£- 1) 



(A.21) 



ALOS'g^ 



2CfI(1 
X — 1] z^ 



99 



X- 1 



{x - 1)2 g4 



zJ Q 



(A.22) 
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(v) e + p^ e + A + X: 

[foDoa,]= Y.e%{x)^D,{z)^oLd',, + Y.elh{x)^D,{z)^oLdi 
q,Q q,g 

+ J2e%{x)AD,{z)AoLa',, , (A.23) 



where 



9,9 
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